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CORRELATION METHOD FOR CHANNEL ESTIMATION FOR OFDM 

The present invention concerns a method and a device for channel 
estimation in communications systems using Orthogonal Frequency-Division 
5 Multiplexing (OFDM). 

The invention is particularly useful in the field of wireless digital 
communications systems. In these systems, one usually deals with multipath 
channels: due to reflections on obstacles, the transmitted signal reaches the 
receiver after having followed a plurality of paths. Hence the received signal 

10 results from the superposition of a plurality of replica of the transmitted signal, 
each one associated (in the equivalent complex baseband representation) with 
a specific delay and attenuation. This received signal is therefore equal to the 
convolution of the transmitted signal with a so-called Channel Impulse 
Response (CIR) /?(f) (where t denotes time). 

"15 The CIR is non-zero only in an interval 0 < t < x, where x is the so- 

called excess delay of the channel, viz., the maximum possible delay between 
the various paths. The channel is said to be time-dispersive if the shape of the 
CIR is different from one narrow peak. Typical values of x range from 50 to 150 
ns for in-door communications, and from 250 ns to 30 \is for mobile radio 

20 communications. 

Equivalents, implementing a Fourier transform of the CIR, the 
time-dispersion properties of a channel may be represented in the "frequency 
domain" by means of a "channel transfer function" H(f); at any given frequency, 
the value of this function is called the "Channel Transfer Factor" (CTF). 

25 The data symbols to be transmitted are modulated, and the receiver 

demodulates the corresponding received signal in order to recover those 
symbols. Clearly, knowledge of the CIR allows one to conduct this 
demodulation more accurately. 

Furthermore, in many modern digital communications systems, for 

30 example telephone systems or mobile radio systems, the demodulation process 
is "coherent". This requires that the system provides a means for "channel 
estimation", viz., for computing an estimate of some channel characteristics of 
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interest (such as the channel impulse response) using a reference portion of 
the received signal. For example, in so-called pilot-assisted channel estimation, 
one transmits a "pilot signal" based on symbols which are a priori known by the 
receiver. As another example, in so-called decision-directed channel 
5 estimation, the receiver uses as a reference some received symbols whose 
value has been determined by means of a tentative decision process. 

A realistic channel always suffers from random noise, which adds up 
to the ideally-noiseless received signal. Therefore, if channel estimation is 
purely based on reference symbols, one only obtains a coarse channel 

1 0 estimate, so that the coherent demodulation process based on this estimate will 
offer sub-optimal performance. In order to improve over such coarse estimates, 
it is often necessary to take into account known statistical properties of the 
channel and/or the noise process. Channel variations over time or over 
frequency have typically a smaller bandwidth than the noise process, so that, 

15 for example, some known channel estimation methods use a filter in order to 
reduce the amount of noise afflicting the channel estimates. 

The need for channel estimation arises, for example, in 
communications systems using Orthogonal Frequency-Division Multiplexing 
(OFDM). 

20 OFDM is a multiple-channel modulation scheme. It is especially 

appropriate for highly frequency-selective channels such as typical channels for 
mobile communications, or for high-rate wireline transmission over copper lines. 
Such channels are characterized by impulse responses which are substantially 
longer than one sample interval. This means that, even in the noiseless case, 

25 each received sample in the digital baseband domain is a superposition of 
multiple transmitted samples weighted by the appropriate channel coefficients. 
In order to resolve such "intersample interference", some kind of equalization 
needs to be performed. 

The way OFDM combats intersample interference is by dividing the 

30 total channel bandwidth into a number D of substantially smaller portions, 
called subchannels. One OFDM channel comprises one parallel use of all 
subchannels. The data to be transmitted is collected into so-called OFDM 
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symbols, and each OFDM symbol is transmitted in parallel on a number D u 
(0 < D u <; D) of these subchannels. The transmitted subchannel signals are 
orthogonal to each other. Since the duration of one OFDM symbol is much 
longer than the sample interval, the problem of intersymbol interference is 
5 strongly reduced. 

To get totally rid of intersymbol interference, a guard interval 
between two symbols is usually introduced during OFDM transmission. If the 
length of the guard interval exceeds the length of the channel impulse 
response, there is no residual intersymbol interference. Furthermore, if the 
guard interval is used in the form of a cyclic prefix, as is usually the case for 
OFDM transmission, one may implement a very simple equalization of the 
frequency-selective channel in the frequency domain. 

The Wireless Local Area Networks (WLAN) systems are examples, of 
radio communications systems which use OFDM. ETSI (European 
15 Telecommunication Standard) BRAN (Broadband Radio Access Network) 
includes a short-range high-data-rate communications system called 
"HIPERLAN type 2" (HIPERLAN/2). HIPERLAN/2 may be used to transport 
Internet Protocol (IP) packets, and will also be capable to act as a wireless 
Asynchronous Transfer Mode (ATM) system, as well as a public access 
20 system, e.g. with an interface to the Universal Mobile Telecommunications 
System (UMTS). The physical layer of HIPERLAN/2 is based on OFDM, with a 
guard interval in the form of a cyclic prefix. Other WLAN systems based on 
OFDM have been standardized by ARIB in Japan (MMAC and its extensions), 
and by IEEE in the US (IEEE802. 1 1 a and its extensions). 
25 ln Pilot-assisted channel estimation schemes for OFDM, known 

symbols are transmitted on given subchannels and time instants for training 
purposes; for example, in systems according to the HIPERLAN/2 or to the 
IEEE802.11a standard, there are two full OFDM pilot symbols preceding every 
burst of data-carrying OFDM symbols. The principle of decision-directed 
30 approaches for OFDM is quite similar: before channel estimation, some data 
symbols are being decided; these decided symbols are then treated in the 
same way as pilot symbols. 
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The simplest channel estimation method consists in comparing a 
transmitted pilot symbol with the received value in the respective subchannel. 
The ratio between these two quantities then yields an estimated subchannel 
transfer factor. This method is known as the "least-squares estimation". 
5 rt ma y be useful to exploit correlations inside the channel transfer 

function, either in the frequency or in the time domain, or both. Assuming, for 
simplicity, that only correlations in the frequency domain shall be exploited in 
the OFDM signal, the characteristics of the channel transfer function over 
frequency can be regarded as a "band-limited" process. Band-limited means 
10 here that the transformation of the channel transfer function into the time- 
domain yields a CIR of limited length, i.e. the CIR length is substantially smaller 
than the OFDM symbol length. This usually holds for OFDM transmission 
systems. 

Therefore, correlations in the frequency-domain may be exploited by 
15 means of an appropriate filter. Ideally, this filter, which is called a "channel 
estimation filter", should be perfectly adapted to the band-limited channel 
frequency response. This means that the filter should be designed in such a 
way that the spectrum of its impulse response displays the same band- 
limitation as the channel frequency response. Hence, an appropriate filter 
20 design requires knowledge of the channel excess delay, which specifies said 
band limitation. 

An example of the use of channel estimation filters to improve the 
quality of channel estimation can be found in the paper by V. Mignone and A. 
Morello titled "CD3-OFDM: A Novel Demodulation Scheme for Fixed and 

25 Mobile Receivers", IEEE Trans. Comm., vol. 44, No. 9, pp. 1144-1151, Sep. 
1996. According to this scheme, which uses a decision-directed coherent 
demodulation, received signals are fed back after they have been corrected by 
means of a channel coding scheme (which needs to be fairly powerful in order 
to ensure stability). Two types of appropriate filters are discussed in this paper. 

30 The first one has a pre-set bandwidth which is flat in the time-domain (and is of 
the order of the guard interval): such a filter is therefore not matched to the 
actual channel excess delay as it optimally should. The second one has a 
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bandwidth which is adaptively determined after the receiver has calculated the 
CIR by Fourier-transforming the CTF: but this conversion from the frequency- 
domain to the time-domain is computationally expensive. 

Besides the purpose of channel estimation filtering, the channel 
5 correlations in the frequency-domain may also be exploited for link-adaptation 
purposes. "Link-adaptation" means that one or more channel-characterizing 
parameters are used to select an appropriate modulation and coding scheme 
for transmission. In particular, the degree of such correlations may serve as one 
such channel-characterizing parameter. 

An example of this approach can be found in the paper by S. Muneta 
et at. titled "A New Frequency-Domain Link Adaptation Scheme for Broadband 
OFDM Systems", belonging to the Proceedings of the 50 th IEEE Vehicular 
Technology Conference (VTC 1999-Fall), pp. 253-257, Amsterdam, The 
Netherlands, Sep. 1999. In this paper, a link-adaptation algorithm is proposed, 
15 in which two pieces of link information are used to adapt coding rate and 
modulation scheme of the current radio link in a HIPERLAN/2-type system. 
These pieces of information are two functions, the values of which depend on 
both the excess delay and the Signal-to-Noise Ratio (SNR) of the channel. 
These functions are calculated from the amplitudes of frequency-domain 
20 samples of the received preamble. But these functions are defined in a more or 
less heuristic fashion, and do not take into account the phase of said 
frequency-domain samples, so that this approach does not perform reliably. 

Thus, according to a first aspect, the invention concerns a method 
for estimating the time-dispersion of a channel comprising D subchannels, 
25 wherein one computes from a received signal a set of estimated Channel 
Transfer Factors (CTF's) H[x>], where v (0<u < D) is the subchannel number, 
said method comprising a step of calculating, for a predetermined strictly 
positive integer d, a correlation factor C d representing the correlations, both in 
amplitude and in phase, between pairs £[o] and H[» + d] of said computed 
30 CTF estimates. 

One should note that the CTF's involved here can be coarse CTF 
estimates, for example least-squares estimates, but they can of course just as 
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well be filtered CTF estimates. Thus, the method according to the invention 
may be adapted to the quality or accuracy of the channel estimation actually 
used. 

Furthermore, the value of parameter d may be adjusted as a function 
5 of the expected order of magnitude of the excess delay x. When the excess 
delay is relatively large, the channel transfer function varies rapidly with 
frequency, and it is therefore necessary to use a small value of d (the minimum 
being 1) in order to keep up with these variations; the averaging over pairs of 
CTF estimates will then include a large number of terms. When the excess 
delay is relatively small, the channel transfer function varies less rapidly and 
one may, therefore, choose larger values of d (2, or bigger) in order to reduce 
complexity, and also in order to increase the quality of the estimation: indeed, 
since a small excess delay leads to a small amount of correlations between 
adjacent subchannels, choosing a bigger d allows one to obtain larger 
15 correlation values and hence to distinguish more easily between different 
excess delays. 

Advantageously, the method according to the invention requires a 
small amount of computation, as can be appreciated on the basis of the two 
following exemplary embodiments. 
20 In one embodiment, said correlation factor C d , when normalized, is 

calculated according to the expression: 
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where the sums over v are carried over available pairs of said computed CTF 
estimates. 

By "available pairs", we mean that pairs where both h[u] and 

H[v + d] have been computed in a prior channel estimation step, either with, or 
without using filters. 
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Whatever the method used to compute the estimated CTF's, they 
will embody a certain amount of noise. This noise has therefore an impact over 
the value of the denominator in the above-stated expression for C d . More 
specifically, this denominator is proportional to the mean received power over 
5 the subchannels belonging to the above-mentioned set of available pairs. This 
mean received power is the sum of a useful received power and of a noise 
power, at least when the noise components are uncorrelated. Hence, this 
denominator depends on the power of the channel estimation error. 

To take into account this effect of the noise, in another embodiment, 
10 said correlation factor C d , when normalized, is calculated according to the 
expression: 



2 

1 



Z#*[»]H[» + d] 
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where C, u is the "mean channel estimation signal-to-noise ratio", and the sums 
over v are carried over available pairs of said computed CTF estimates. 
15 This mean channel estimation signal-to-noise ratio <; u , which 

specifies the quality of the channel estimates, is defined as: 



<,u =— - , 



where: 

P h =s{\H[»]\ 2 \ 

20 is the mean-squared channel amplitude per subchannel, and 
Pe=e§H[v]-H[»]\ 2 \ 

is the mean power of the effective subchannel estimation error. 

The method according to the invention can be useful for many 
applications where knowing the time-dispersion characteristics of a channel is 
25 required, and is, for example, particularly suitable for designing a channel 
estimation filter, and for link adaptation. 
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In one such application, the method further comprises a step of 
looking-up in a pre-constructed mapping table a value of channel excess delay 
x corresponding to the value of said correlation factor C d . 

This mapping table is constructed by associating each value of C d 
5 displayed in it with the corresponding value of channel excess delay measured, 
in one way or another, in the same physical setting as the one in which this 
value of C d has been calculated. Equivalently, one may use an experimentally 
established semi-analytic relation instead of a mapping table. 

In another application, the method further comprises a step of 
10 adapting some link parameters as a function of the value of said correlation 
factor C d . These link parameters may include for example the modulation 
scheme and/or the coding scheme. 

According to a second aspect, the invention concerns a device for 
estimating the time-dispersion of a channel comprising D subchannels, said 
15 device receiving as an input a set of estimated Channel Transfer Factors 
(CTPs) H[v], where v (0 ^ u < D) is the subchannel number, computed from a 
received signal, characterized in that it comprises a correlations unit capable of 
computing a correlation factor C d , where of is a predetermined strictly positive 
integer, representing the correlations, both in amplitude and in phase, between 
20 pairs H[»] and H[o + d] of said computed CTF estimates. 

In a particular embodiment, this time-dispersion estimation device 
also comprises a parallel-to-serial unit capable, when provided with a CTF 
vector H as an input, of providing said correlations unit with a series of 
individual CTPs h[o] classified by successive subchannel number v. 

25 This device may use various appropriate expressions for said 

correlation factor C d . 

In one embodiment, In one embodiment, said correlation factor C d , 
when normalized, is calculated according to the expression: 
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where the sums over v are carried over available pairs of said computed CTF 
estimates. 

In another embodiment, where the effect of the noise is taken into 
5 account, the correlation factor C d , when normalized, is calculated according to 
the expression: 



where <; u is the mean channel estimation signal-to-noise ratio, and the sums 
over v are carried over available pairs of said computed CTF estimates. 

10 The time-dispersion device according to the invention may 

conveniently be applied to the estimation of the excess delay of a channel. For 
this purpose, this device will also comprise a look-up table, capable of providing 
a value of channel excess delay x corresponding to the value of C d . 

The time-dispersion device according to the invention may also 

15 conveniently be applied to the estimation of link quality. For this purpose, this 
device will also comprise a link adapter responsive to the value of said 
correlation factor Cd. 

The invention also relates to: 

- a modulated-signal reception apparatus including a device as 
20 succinctly described above, 

- a telecommunications network including at least one reception 
apparatus as succinctly described above, 

- a data storage means containing computer program code instructions 
for executing steps of anyone method as succinctly described above, and 

25 " a computer program containing instructions such that, when said 

program controls a programmable data processing device, said instructions 
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mean that said data processing device implements anyone method as 
succinctly described above. 

The advantages offered by these devices, apparatus, 
telecommunications networks, data storage means and computer programs are 
5 essentially the same as those offered by the methods according to the 
invention. 

Other aspects and advantages of the invention will emerge from a 
reading of the detailed description, which will be found below, of preferred 
embodiments given by way of non-limitative example. This description refers to 
10 the accompanying drawings, in which: 

- Figure 1 is a diagram illustrating the standard time-domain structure of 
an OFDM preamble, and 

- Figure 2 depicts schematically a time-dispersion estimation device 
according to the invention. 

15 For introduction, consider the transmission of one single OFDM 

symbol generated with a D-dimensional Inverse Discrete Fourier Transform 
(IDFT). Each of the D subchannels is modulated by a complex-valued 
amplitude A M [u], where v is the subchannel number (0Zv<D) representing 
subchannel frequency, and u indicates the OFDM symbol number representing 

20 transmission time. The data are mapped on the subchannel amplitudes A ]p]. 
The OFDM symbol is equipped with a guard interval of D g samples. The 
discrete-time complex baseband OFDM transmitted signal a M for OFDM symbol 
number \x, when written out as a vector over discretised time p, with: 
- D g <> p < D, has the following components: 

which can be generated by an IDFT of dimension D. 

Using s to denote expectation value (computed over a set of 
samples), the average signal power is 
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where E s is the average energy per channel symbol and T is the modulation 
interval. 

For convenient notation, the symbol index n will be suppressed in the 

following. 

5 We assume transmission over a multipath channel, whose CIR will 

be modeled, during the transmission of each OFDM symbol, by means of a 

discrete-time vector h having components h[ P ] (p = 0 D-1). The fact that the 

continuous-time CIR is non-zero only up to the excess delay x is equivalent to 
having only the components h[ P ] with p = 0,...,D e -1 (where D e is the discretised 
1 0 excess delay of the CIR) being non-zero. 

As discussed in the introduction, the noiseless received signal is 

given by 

y=a*h, 
where * denotes circular convolution. 
1 5 For convenient notation, without loss of generality, perfect frequency 

synchronization will be assumed in the receiver. At the receiver input, samples 
n[p] of Additive White Gaussian Noise (AWGN) corrupt the noiseless signal 
samples tfp] to yield the actually-received samples 
y[p] = ftp] + n[p]. 
20 The noise power is 

where N 0 is the "one-sided power spectral density" of the white noise. The 
channel Signal-to-Noise Ratio (SNR) at the receiver input is thus 

Implementing a Discrete Fourier Transform (DFT), one defines the 
discrete Channel Transfer Factor (CTF) at subchannel u as 



25 
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The OFDM receiver yields the noisy subchannel amplitudes via DFT 
processing as 
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D-l ,2*. 

E:" " 

= f M + Nfo] , 

where 

— A[v] h[»] (3) 
5 is the noiseless received subchannel amplitude, and N[»] is the resulting noise 
sample on subchannel u. Here we assumed that the guard interval D g is 
sufficiently large so that the linear convolution of the transmitted OFDM signal 
and the CIR is equal to their circular convolution in the time interval of interest. 

It is obvious from Equation (3) that for data recovery from y[o], an 
1 0 estimate of h[u] is required. 

We shall now describe, as an example of a method which can be 
implemented in the receiver to obtain a coarse estimate of H[v] within the 

framework of the present invention, the known "least-squares" estimation 
method. 

1 5 The active subchannel positions will be referred to by the subchannel 

indices v which belong to a set U containing the D u indices of these positions. 
The remaining (D - D u ) subchannels, namely the ones in which the modulation 
amplitude is set to zero, are usually called "virtual" subchannels. 

We shall assume that, for so-called training purposes, a preamble 

20 symbol containing known subchannel amplitudes covering the entire active 
frequency-domain is provided; we will collect these amplitudes in a vector A p of 
length D u . 

The transmitted time-domain preamble (or pilot) OFDM symbol is 
depicted in the upper part of Figure 1. The repeatedly transmitted time-domain 
25 sequence of length D is obtained by the IDFT 

1 

ap = VF Dl,tAp ' W 

where we introduced the D u xD partial DFT matrix 
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and f denotes the conjugate transpose of a matrix. 

The sample sequence resulting from the convolution of the 
transmitted preamble sample sequence with the CIR h is disturbed by additive 
5 noise. Since the guard interval and the repeated preamble, which are assumed 
to be of sufficient size, are both cyclic, the linear convolution is equal to the 
circular convolution for these signal segments which are of interest for final 
demodulation. Thus, the noiseless received sequence y p , of length D is equal to 
the circular convolution of a p and h. 

10 ,n tne ,ower Part of Figure 1, the temporal structure of the received 

signal is indicated. Perfect channel frequency synchronization has been 
assumed. The successive received sequences y u y 2 , and so on, are made of 
the repeated noiseless sequence y p disturbed respectively by the (statistically 
independent) successive D-dimensional noise vectors n 1f n 2 , and so on. 

15 Implementing a DFT, we find: 

1 1 
Y P s ^ D «yp =d ' a g{^ D u a p }- D u h=diag{Ap}H = AH , (6) 

where we collected the D u channel transfer factors H[u] (v e U) in the column 
vector 

H = D u h , ^ 

20 and we collected the transmitted pilot subchannel amplitudes in the D u x D u 
diagonal matrix 

A = diag { Ap } . (8) 

Each received sequence is of the form: 

yp = Yp + n . (9) 
25 where n is a noise vector. 

This received signal y p is transformed into frequency domain via the 



DFT 



Y P s ^ D uyp = A-H+_ / L D u n = A H + N , ( 10 ) 
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where we used Equation (6), and we introduced the D u -dimensional frequency- 
domain noise vector N. 

Finally, the "least-squares estimate" H LS of the CTF vector is 

obtained by dividing out the frequency-domain received signal by the pilot 
5 subchannel amplitudes: 

h ls=A" 1 Y p . ^ 

This operation can be seen as extracting or removing the influence 
of the transmitted pilot amplitudes from the received signal. 

The least-squares estimation method just described is only one 
10 possible way to obtain a coarse estimate H coarse of the CTF. Furthermore, for 
good quality channel estimation, filtering of H coarse by means of a channel 
estimation filter is commonly arranged in the receiver. But the present invention 
teaches a method which allows one to estimate the channel excess delay x 
whether the measured channel transfer factors have, or have not, been filtered 
15 beforehand; whatever the case may be, these CTF estimates will be denoted 
as H\o) in the following. 

A characteristic to quantify the amount of time-dispersion will now be 
derived. This characteristic exploits the close connection existing between the 
size of the excess delay, representing essentially the width of the CIR, and the 

20 amount of self-correlations in the channel transfer function, since the latter is 
the Fourier transform of the CIR (as expressed by Equation (2) above). 

This auto-correlation is explored, according to the invention, by 
comparing, both in amplitude and in phase, the channel transfer factors 
corresponding to subchannels separated by a predetermined subchannel 

25 number o*. The amplitude of the difference between such pairs is summed over 

"available pairs", viz., pairs where both H[»] and H[» + d] have been computed 

by some channel estimation means, and preferably over all such pairs; this sum 
shall be written as: 

s["hT\Hb + d]-e +J *rt»]\ 2 . (12) 
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where the constant <p, the so-called "phase ramp", takes into account a possible 
shift in the time window determining, for each OFDM symbol, the set of 
samples involved in the DFT (see Equation (10) above). 

The phase ramp contribution to s|hJ in Equation (12) is of no interest 
as far as channel dispersion estimation is concerned. One can get rid of this 
phase ramp, and hence of the dependence upon potentially different time 
windows, by minimizing the expression of s[h] with respect to <p: 

( r 
min s[H\ = min £|£[o + d] | 2 + £|£[o] | 2 _ 2 . Re h.-»A*fp]A\p + d] 



J 



(13) 
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-2|^+4| 2 +Z|A[o]| 2 -2. 

If the CIR exhibits, a large dispersion,.-then the frequency response 
will be very selective, so that this difference power will be large. Obviously, the 
characteristic depends on the absolute power of the CIR. It is therefore 
reasonable to factor out this CIR power, to obtain: 



min s[h\ =h \H[„ + d ]\ 2 + £|£[o] I 2 1 (1 _ Q ) , 

V lu u J 

1 5 where we introduced a normalized correlation factor 



(14) 



C d s 



Z^[o]| 2 + |^[oW]| 2 ) 



(15) 



20 



A value of C d near one indicates a rather non-selective CIR with 
small channel power dispersion, whereas smaller values of C d characterize 
larger dispersion. 

This dispersion estimation can be improved upon by taking into 
account the fact that the estimated CTF's may embody a certain amount of 
noise. As explained above, this noise has an impact over the denominator 
(which plays the role of a normalization factor) in Equation (15). It is however 
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possible to get rid of this noise effect on the normalization factor, as will be 
shown now. 

Let us consider the most simple case of a flat channel, and of real- 
valued subchannel transmitted symbols, e.g. a so-called BPSK modulation. 
5 Then, the expectation value of C d is 

p H _ 1 



P H+P* 1 + J_ 



where P H , p e , and £ u have been defined in the introduction above. Since in the 
case of a flat channel, one expects a correlation factor equal to 1, it appears 
that one should use, in this case and hence in general, 



10 C d = 



\ + ±] 



^H*[v]H[v + d] 



x{\m\ 2 +\a[»+d]\ 2 ) 



(16) 



as a noise-corrected normalized correlation factor. 

One should note that, in the above derivation, the subchannel 
estimation errors f #[v]-#[o]|} were assumed to be uncorrelated. Therefore, 

the mean channel estimation signal-to-noise ratio <; u must be evaluated before 
15 any frequency-domain filtering, since such filtering produces correlated noise. 

For example, in the case the estimation is carried out by means of 
the least-squares estimation method as described above, the mean channel 
estimation signal-to-noise ratio is equal to: 

^u=(D/D u K. 

20 This expression reflects the fact that the signal power is concentrated in the 
active subchannels. 

On the other hand, time-domain filtering (for example implemented 
by averaging over two successive preamble sequences yi and y 2 in the time- 
domain) does not cause any frequency-domain noise correlation. When one 

25 uses such a filtering, it should be taken into account when evaluating <;„. One 
gets in this case: 
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C u s2(0/D U K, 

where the factor of 2 reflects the noise reduction obtained thanks to the 
averaging process. 

Figure 2 depicts schematically an exemplary time-dispersion 
5 estimation device 100 according to the invention. 

The input to this device 100 is a set of CTF's which have been 
computed from the received signal, possibly collected in the form of a CTF 
vector H. This set is processed by a parallel-to-serial unit 101, whose output 
consists of the series of individual CTF's (the components of H) classified by 
10 successive subchannel number. 

Next, correlations unit 102 computes the correlation factor C d 
according to Equation (15). Optionally, the correlation factor C d may be 
computed according to Equation (16), in which case the value of the mean 
channel estimation signal-to-noise ratio ^ u is also provided to correlations unit 
15 102. This unit comprises delay registers adjusted to the value chosen for the 
integer parameter d, so as to be able to associate in the calculation pairs of 
CTF's whose subchannel numbers differ by d. 

Finally, correlations unit 102 transmits the calculated value of the 
correlation factor C d to a Look-Up Table (LUT) 103, in which the value of 
20 excess delay x corresponding to this value of C d is read or interpolated, before 
being sent to the output of the time-dispersion estimation device 100. 

As a variant, look-up table 103 may be replaced by a calculation unit 
capable of calculating x when given C d by means of some experimentally 
established semi-analytic relation. 



